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We present a multipartite entanglement measure for iV-qubit pure states, using the 
norm of the correlation tensor which occurs in the Bloch representation of the state. We 
compute this measure for several important classes of iV-qubit pure states such as GHZ 
states, W states and their superpositions. We compute this measure for interesting appli- 
cations like one dimensional Heisenberg antiferromagnet. We use this measure to follow 
the entanglement dynamics of Grover's algorithm. We prove that this measure possesses 
almost all the properties expected of a good entanglement measure, including monotonic- 
ity Finally, we extend this measure to iV-qubit mixed states via convex roof construction 
and establish its various properties, including its monotonicity. We also introduce a re- 
lated measure which has all properties of the above measure and is also additive. 

PACS numbers: 03.67.Mn, 03.65. Ca, 03.65.Ud 



I. INTRODUCTION 



Entanglement has proved to be a vital physical resource for various kinds of quantum 
information processing, including quantum state teleportation [1,2], cryptographic key 
distribution [3], classical communication over quantum channels [4,5,6], quantum error 
correction [7], quantum computational speedups [8] and distributed computation [9,10]. 
Further, entanglement is expected to play a crucial role in the many particle phenomena 
such as quantum phase transitions, transfer of information across a spin chain [11,12] etc. 
Therefore, quantification of entanglement of multipartite quantum states is fundamental 
to the whole field of quantum information and in general, to the physics of multicomponent 
quantum systems. Whereas the entanglement in pure bipartite states is well understood, 
classification of multipartite pure states and mixed states, according to the degree and 
character of their entanglement is still a matter of intense research [13,14.15]. Principal 
achievements are in the setting of bipartite systems. Among these, one highlights Woot- 
ter's formula for the entanglement of formation of two qubit mixed states [16], which still 
awaits a viable generalization to multiqubit case. Others include corresponding results 
for highly symmetric states [17,18,19]. The issue of entanglement in multipartite states is 
far more complex. Notable achievements in this area include applications of the relative 
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entropy [20], negativity [21] Schimidt measure [22] and the global entanglement measure 
proposed by Meyer and Wallach [23] . 

A measure of entanglement is a function on the space of states of a multipartite system, 
which is invariant on individual parts. Thus a complete characterization of entanglement is 
the characterization of all such functions. Under the most general local operations assisted 
by classical communication (LOCC), entanglement is expected to decrease. A measure 
of entanglement which decreases under LOCC is called an entanglement monotone. On 
bipartite pure states the sums of the k smallest eigenvalues of the reduced density matrix 
are entanglement monotones. However, the number of independent invariants (i.e. the 
entanglement measures) increase exponentially as the number of particles A" increases and 
complete characterization rapidly becomes impractical. A pragmatic approach would be 
to seek a measure which is defined for any number of particles (scalable), which is easily 
calculated and which provides physically relevant information or equivalently, which passes 
the tests expected of a good entanglement measure [13,14]. 

In this paper, we present a global entanglement measure for A^-qubit pure states which 
is scalable, which passes most of the tests expected of a good measure and whose value 
for a given system can be determined experimentally, without having a detailed prior 
knowledge of the state of the system. The measure is based on the Bloch representation 
of multipartite quantum states [24] . 

The paper is organized as follows. In section II we give the Bloch representation of 
a A-qubit quantum state and define our measure Er- In section III we compute Er for 
different classes of V-qubit states, namely, GHZ and W states and their superpositions. 
In section IV we prove various properties of Er, including its monotonicity, expected of 
a good entanglement measure. In section V we extend Er to V-qubit mixed states via 
convex roof and establish its monotonicity. Finally, we conclude in section VI. 

II. BLOCK REPRESENTATION OF A V-QUBIT STATE AND THE 
DEFINITION OF THE MEASURE 

Consider the generators {I,o~ x , a y , a z } = {oo, o~2, er 3 } of SU(2) group (Pauli ma- 
trices). These hermitian operators form a orthogonal basis (under the Hilbert-Schmidt 
scalar product) of the Hilbert space of operators acting on a single qubit state space. The 
N times tensor product of this basis with itself generates a product basis of the Hilbert 
space of operators acting on the A^-qubit state space. Any A-qubit density operator p can 
be expanded in this basis. The corresponding expansion is called the Bloch representation 
of p [24]. 

In order to give the Bloch representation of a density operator acting on the Hilbert 
space C 2 <g)C 2 ®- • -<g>C 2 of a A^-qubit quantum system, we introduce following notation. We 
use k, ki (i — 1, 2, • • • ) to denote a qubit chosen from A" qubits, so that k, ki (i — 1, 2, • • • ) 
take values in the set M = {1,2,- •• , A^}. The variables or for a given k or 
ki span the set of generators of SU{2) group for the kth or kith qubit, namely the set 
{I in, o\ k , <72 k . , o~3 k . } for the /c;th qubit. For two qubit k\ and A>2 we define 
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= (J 2 <g> J 2 <8> • • • ® a aki <g> I 2 <g> • • • <g> h) 

<T%£ = (I 2 <g> h ® ■ • • ® ^ 2 ® / 2 ® • • • ® / 2 ) 

= ( J 2 ® / 2 ® • • • ® <T afcl <g> J 2 ® ■ ■ ■ ® V ak2 ® h ® / 2 ) (1) 

where cr afcl and <7 Qfc2 occur at the fcith and fc 2 th places (corresponding to fcith and 
/c 2 th qubits respectively) in the tensor product and are the a^th and «fc 2 th generators of 
£77(2), (a kl = 1, 2, 3 and a k2 = 1, 2, 3) respectively. Then we can write 

p = ^fafr + E E w£> + E E + • • ■ + 

{k}cAf "fe {fci,fc 2 } Qfe 1 afc 2 

V V t a (ki) a (k2) (k M ) + .,. + \^ t (1) (2).. (JV)| 

{fci,fc 2 ,- ,&Af} a fc 1 Q fe 2 '" a feM aia 2 ---ajv 

(2) 

where s^) is a Bloch vector (see below) corresponding to fcth subsystem, s( fe ) = 
[s Qfc ]a fc =i which is a tensor of order one defined by 

s ak = Tr[pa^} = Tr[p k a ak ], (3) 

where pk is the reduced density matrix for the fcth qubit. Here {ki,k 2 ,--- , /cm}, 1 < 
M < N, is a subset of Af and can be chosen in (^) ways, contributing (^) terms in the 
sum X^{fcx fc 2 - fc M } * n Eq-(2), eacn containing a tensor of order M. The total number of 
terms in the Bloch representation of p is 2^. We denote the tensors occurring in the sum 
E {kuk2 ,.., kMV (1 < M < N) by T^-M = [t akiak2 ... ak J which are defined by 

l<a kl a k2 ...a kM \P U a kl u a k2 a k M ' 

= Tr\p klk2 ... kM (a aki <g> a ak2 <g> • • • <g> <r QfcM )] (4) 

where p klk2 ... kM is the reduced density matrix for the subsystem {fci/cv- • • • k M }. We call 
The tensor in last term in Eq. (2) T^ N \ 

From Eq.(4) we see that all the correlations between M out of N qubits are contained 
in 7~{ fe i> fe 2,---,fcM} anc } a n the _/y qubit correlations are contained in T^ N \ 

If p is a iV-qubit pure state we have 

1 N 

Trp 2 = ^r(l + ^ ||s (k) || 2 + ^ ||T^- fc2 }|| 2 + ---+ ^ || r {*i,fc2,-,fcM}||2 

k=l {fci,fe 2 } {fei,fe2,-- >^m} 

+ ... + ||TW|| 2 ) = l (5) 
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Any state p = \ip){ip\ living in d 2 dimensional Hilbert space of operators acting on 
a d dimensional Hilbert space of kets, can be expanded in the basis comprising d 2 — 1 
generators of SU(d) and the identity operator. The set of coefficients in this expansion, 
namely {Tr(pXi)} i = 1, 2, • • • , d 2 — 1 is a vector in M. d ~ l and is the Bloch vector of p. 
The set of Bloch vectors and the set of density operators are in one to one correspondence 
with each other. The set of Bloch vectors for a given system forms a subspace of M d2_1 
denoted B(R d The specification of this subspace for d > 3 is an open problem [25,26]. 
However, for pure states, following results are known [27]. 



II* = «» 

D r (R d2 ~ 1 ) C B(R d2 ~ 1 ) C D^W* 2 - 1 ) 



where D r and Dr are the balls of radii r = \J 2 (d~i) anc ^ ^ = \J ^2 ^ respectively in 
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We propose the following measure for a iV-qubit pure state entanglement 

£t(I^» = (II^ (a °II-i) (7) 

where is given by Eq.(4) for (M = N) in Bloch representation of p = \ip)(ip\- The 
norm of the tensor appearing in definition (7) is the Hilbert-Schmidt (Euclidean) 
norm ||T^|| 2 = (T^ N \T^) = ^2 aia2 ... aN t 2 tlCt2 ... aN . Throughout this paper, by norm, 
we mean the Hilbert-Schmidt (Euclidean) norm. We comment on the normalization of 
E T {\^)) below. 



III. GHZ AND W STATES 

Before proving various properties of Er(\ip)), we evaluate it for states in the iV-qubit 
GHZ or W class. A general A-qubit GHZ state is given by 

|V) = Vp|000---0) + v/l-p|lll---l); A>2 (8) 

A general element of T (JV) is given by t hi2 ... iN = (ip\a h <g> a h <g> • ■ ■ ® cn N \i p), %k = 
1,2,3; k = 1,2, • • • , N. The nonzero elements of T (JV) are t U -i = 2y/p(l-p), t 33 ... 3 = 
p + (—l) N (l—p). Other nonzero elements of are those with 2k a 2 s and (N — 2k) a±s, 
k = 0, 1, • • • , LyJ where [^J is the greatest integer less than or equal to x. (eg. for 
N = 3, t 12 2 etc). These are equal to (— l) k 2^p(l — p). This gives 

Lf J 

\\ T(N) 1 1 2 = Ml - P ) + (p+ (-1)^(1 - P)) 2 + 4p(l -p)J2 (1) ■ (9) 



4 



Thus we get, for Er{\ip)) 




Eq.(8), with N = 2, represents a general two qubit entangled state in its Schmidt 
decomposition 

i^) = vp|oo> + v/r^iii). 

Thus Eq.(lO) gives the entanglement in a two qubit pure state. Using Eq.(lO) it is 
straightforward to see that Eq-(\ip)) for an arbitrary two qubit pure state is related to 
concurrence by 

£r(M) = Vl + 2C 2 -1, 

where concurrence C for such a state is 2^p(l — p). 

Figure 1 plots E T (^)) in Eq. (10) as a function of p for N = 3. For the iV-qubit GHZ 
(maximally entangled) state p = 1/2, so that 




We see that, as a function of N, Eq-(\GHZ)) increases as a polynomial of degree 
LfJ. Figure 2 plots E T (\GHZ)) as a function of N. E?(\GH Z)) increases sharply 
with N as expected. Note that Er(\GH Z)) > for GHZ class of states. Whenever 
appropriate, we normalize the entanglement of a iV-qubit state Er(\ip)) by dividing 
by R N = E T {\GHZ)). 

iV-qubit W state is given by 




00---L,0---00); N>3 



where jth term has a single 1 at the jth bit. The state \W) = ^ =1 a[ \W) is given by 
\W) = ^j- 1 11 • • • Ojl • • ■ 11) ; N > 3 and has a single at the jth bit. We note that 

\W) is locally unitarily connected to \W) so that their entanglements must measure to 
the same value. The general element of for the state p = \ W){W\ is 

1 N 
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FIG. 1. Variation of Eq-QGHZ)) for N = 3, express in units of R3, with parameter p. 
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FIG. 2. Variation of E r {\GHZ)) with number of qubits N. 



+ N (00---l r --00| ( j il ®(7 i2 ®---®(7 ijv |00---l r --00) 
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Only the first term contributes to £33... 33 = —1. Other nonzero elements have the form 

*3-31j3-31i3-3 — Jf — ^3---32 j 3---32 ; 3---3- 

There are (^) elements of each of these two types, so that 

m = 1 + 2 Q)'Q =1 + 4 ^l (12) 

E?(\W)) = WT^W -1 = ^1 + 4^1-1 (13) 

It is straightforward to check that Er(\W)) = Er(\W)) as expected. Note that 
E r {\W)) > 0. 

Next we consider a superposition of \W) and \W) states, \ip s ,<i>) — y/~s\W)+\/l — se l ^\W). 
It is clear that the entanglement of \ip s ,4>) cannot depend on the relative phase 0, as \ip s ,(f>) is 
invariant under the local unitary transformation (|0), |1)} — > {|0), e l<?i |l)} upto an overall 
phase factor. As we shall prove below, Er is invariant under local unitary transforma- 
tions. Figure 3 shows the entanglement of \ip s ,<t>) as a function of s, calculated using our 
measure. 



0.95 - 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

s 

FIG. 3. Variation of Eq-(\ip s ,(p)), expressed in units of R N , with the superposition 

parameter s. 



An important example of W state and its generalizations is the one dimensional spin-| 
Heisenberg antiferromagnet, on the lattice of size N, with periodic boundary conditions, 
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given by the Hamiltonian 



N 

H N = ^2XjX j+1 + YjY j+1 + ZjZ j+1 



(14) 



where the subscripts are mod N and X, Y, Z denote Pauli operators a x , a y , a z respectively. 
Hn commutes with S z = Y^ Zj, so the eigenstate of Hn is a superposition of basis vectors 
|&i--- b n ) with s of b-i ■ ■ ■ b N are ones and N — s are zeros for some fixed < s < N. When 
s = 1, the translational invariance of Hn implies that the eigenstates are 



7V-1 



5^e <fcj "|00---lj0---0> 



(15) 



i=o 



where the jth summand has a single 1 at jth bit just like W state and the wave 
number k = for some integer < m < N — 1. The state is locally unitarily 

transformed to the W state so that it has the same value of E r (\W}) or E r (W)). 
For s > 2 the eigenstates of H N have the form 



l^jv(s)) = 



(J {ii-i4 



^ |oo...i j o-..i Js o), 



(16) 



where 1 occurs at j\ • • ~js, {ji • • • js} Q Af = {1, 2, • • • , N} and can be chosen in (^) 
ways. We see that for ^^(s)), *33-3 = 1- For even TV, ti... 12-23-3 with x 1 s and y 2 s, 
corresponding to the average of x <r x s, y a y s and iV — x — y a z s, we get, for even x and 
even y, 

x\fy\ fx + y\]fN 



*1. 



12-23-3 



x - y 



Since |-0jv(s)) is a symmetric state, any permutation of its indices does not change the 
value of an element of [24], so that, 



2s 

E 

z+y=2 
even 



x + y 



N — x — y 



s — 



N - x 

y 



— 1 with s. We see that it is 



Figure 4 shows the variation of Et(\^n(s))) = II^S S )\I 
maximum at s = y, which is the ground state of H N , as expected. Note that EVG^aKt))) 
for the ground state (s = y) rises far more rapidly than the entanglement of the iV-qubit 
GHZ state _Rjv = Eq-{\GHZ)), (Eq.(ll)), with the number of spins (qubits) N. This can 
be understood by noting that IV'jv(t)) ^ or s = T can ^ e wr itten as a superposition of 
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K— ) ^ _( l u bit GHZ states. For example, |-0 4 (2)) can be written as the superposition of 
three 4-qubit GHZ states, 



1^4(2)) 



1 



L ^=(|0011) + |1100)) + ^(loioi) + |1010)) + ^(liooi) + |0110)) 



As N increases, initially Et(\iI>n(%))) is comparable to Rn, but after N = 16 the ratio 
1 ^- increases very rapidly, reaching 10 7 for 100 qubits. Also, as N increases, 



Ha 



Et(\i/jn(s))) falls off more rapidly as s deviates from y. We are presently trying to 
understand this behavior. 




2 4 6 8 10 12 14 16 18 20 



FIG. 4a 



Finally, in this section, we consider the superpositions of W and GHZ states, 

\^w+ghz(s, 0)) = V~s\GHZ) + VT^~s e^\W), (17) 
also considered in [28]. For three qubits, N = 3, a direct calculation gives, for this state, 



|| T (A0|p 



4s 2 + 6s(l-s) + — (s-1) 2 ; (0<s<l) 
3 

-N 



(18) 



E r (\ip w +GHz(s,(j))}) = \\ r w N ( s ))\\ ~ 1 
which coincides with the corresponding values of W (s = 0) and GHZ (s = 1) states. 
Note that Eq-(\ipw+GHz{s,4>))) is independent of the phase 0, in contrast to the entan- 
glement measure used in [28]. Figure 5 shows the dependence of Er(\ipw+GHz(s, <f>))) on s. 
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10 20 30 40 50 60 70 80 90 1 00 

s 



FIG. 4b 

FIG. 4. Variation of Ej-(\ippf(s))), in units of Rn with s, for (a) N = 20 and (b) 

N = 100 (see text). 



1.02 




s 



FIG. 5. Variation of Er(\ipw+GHz(s, </>))), expressed in units of Rjy, with the 
superposition parameter s, for N = 3. 



10 



IV. PROPERTIES OF E r (\ip)) 

To be a valid entanglement measure, Er{\ip)) must have the following properties [29,30]. 

(a) (i) Positivity : Eq-{\ip)) > for all iV-qubit pure state (ii) Discriminance: 
Er{\i))) = if and only if is separable (product) state. 

(b) LU invariance : Eri\ip)) is invariant under local unitary operations. 

(c) Monotonicity : local operators and classical communication (LOCC) do not in- 
crease the expectation value of E r (\ip)). 

We prove the above properties for Eq-(\ip))- We also prove the following additional 
properties for Et(\i/j)). 

(d) conUnmty - 10X01)11 - \e(\^)) - £(|0»| - 0. 

(e) superadditivity E T (\^) <g> \(f>)) > E r (\^)) + E r (\(f>)). 

We need the following result which we have proved in [24] . 

Proposition : A pure iV-partite quantum state is fully separable (product state) if 
and only if 

TW = sWos( 2 »o...osW, (19) 

where s( fc ) is the Bloch vector of /cth subsystem reduced density matrix. 

The symbol o stands for the outer product of vectors defined as follows. 

Let uW, u®, . . . , u( M ) be vectors in R d ?-\ R^' 1 , ■ ■ ■ , R^m" 1 . The outer product o 
u^ 2 ) o • • • o u( M ) is a tensor of order M, (M-way array), defined by 

l W2-lM — U ii U i 2 ■ ■ ■ U i M ' 1 — lk — ^ J"? K — - 1 ? Z ! J iW " 

Proposition 1 : Let j^) be a iV-qubit pure state. Then, = 1 if and only if 

is a separable (product) state. 

Proof. By proposition 0, \ip) is separable (product) if and only if 

= s^ 1 W 2 ) o ■ ■ • o S W, 

As shown in [31,32 ], 

(OL^ k \OLiS W ) = ^i(s^,s^), (20) 
where (, ) denotes the scaler product. This immediately gives, for qubits, 

I| T (ao||2 = (r w rW) = nf =1 ( s ( fc ), s «) = n fc || s ( fc )|| 2 = i 

Proposition 1 immediately gives 

Proposition 2 Let \ip) a iV-qubit pure state. Then Eq-(\ip)) = if and only if \ip) is 
a product state. 

Proposition 3 : Let be a iV-qubit pure state. Then ||T^|| > I. 
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It is instructive to show this result by direct computation of for the case of two 
and three qubit states. First, consider a general two qubit state 

|V>) = ai|00) + a 2 |01) + a 3 |10) + a 4 |ll), ^ W? = 1 

k 

by direct computation we get 

||T( 2 )|| 2 = l + 8(|a 2 a 3 |-|aia4|) 2 >l (21) 

This means, via proposition 2 that is a product state if |a 2 a 3 | = |aia4|. Next, 
consider a three qubit state in the general Schmidt form [33] 

|^) = A |000) + Aie^|100) + A 2 |101) + A 3 |110) + A 4 |lll) (22) 

where \ > 0, % = 0, 1, 2, 3, 4 and £\ A 2 = 1 
By direct calculation of we get 

||T (3) 1 1 2 > 1 + 12A 2 A 2 + 8A 2 A 2 + 8A 2 A 2 + 8(AxA 4 - A 2 A 3 ) 2 > 1 (23) 

Here the conditions for product state become (AiA 4 = A 2 A 3 and A = 0). We now 
prove proposition 4 for a general iV-qubit state \ip). 

If is not a product of N single qubit states (i.e. is not ^-separable) then it is 
N — k separable k — 2, 3, • • • ,N — 1. Viewing the iV-qubit system as a system comprising 
N — k qubits, each with Hilbert space of of dimension 2 and k entangled qubits with the 
Hilbert space of dimension 2 k , we can apply proposition to this separable system of 
N-k + 1 parts in the state \^). We get 7jg° = (s«) o (s^) o • • • (s iN ^) o (s^-^ 1 )) 

This implies, as in proposition 1, via Eq. (20) and Eq. (6) that 

ll^ll = nf_7 fc+1 ||(^)|| 2 = > i (4 = 2% (24) 

If k = N we attach an ancilla qubit in an arbitrary state \<j>) and apply proposition to 
iV + 1 qubit system in the state <8> \4>) where \ip) is the A^-qubit entangled state. This 
result, combined with proposition 1, completes the proof. 
Proposition 4 immediately gives 

Proposition 5 : E<r{\ip)) > 0. 

We now prove that Eq-(\ip)) is nonincreasing under local operations and classical com- 
munication. Any such local action can be decomposed into four basic kinds of operations 
[34] (i) appending an ancillary system not entangled with the state of original system, 
(ii) performing a unitary transformation, (iii) performing measurements, and (iv) throw- 
ing away, i.e. tracing out, part of the system. It is clear that appending ancilla cannot 
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change ||T^||. We prove that Er(\i>)) does not increase under the remaining three local 
operations. 

Proposition 6 : Let C/j (i — 1, 2, • • • , N) be a local unitary operator acting on the 
Hilbert space of ith subsystem TiS^. 
Let 

P = ^liU^p'^lM) (25) 

for density operators p and p' acting on 7i = ®f =1 7i^ and let 1~( N ) and denote the 

N partite correlation tensors for p and p' respectively. Then, 
|| T /(A0|| = \\TW\l so that E r (p) = E T (p') 

Proof. Let U denote a one qubit unitary operator, then it is straightforward to show 
that Ua a W = Y^pOapo-p 

where [O a p] is a real matrix satisfying 00 T = I = T 0. It is an element of the 
rotation group 0(3). Now consider 

4*2-^ = Tr (p' a h ® ® • • • ® o- iN ) 

= Tr^^Mo-n ® a l2 ■ ■ ■ a iN (®tiUl)) 
= TripU^Ul <g> U 2 a i2 U\ <g> • • • <g> U N a iN U^ N ) 



J2 Tr(pa ai ® a a2 ® • • • ® a Q JO^O^ • • • O} 

= V t o (1) o (2) • ••o (A0 

ai---Qjv 



= (T^x 1 0^x 2 0^---x N 0^) tll2 ... tN 

where is the /c-mode product of a tensor e jj3x3x---3 D y or thogonal matrix 
(fe) e R 3x3 [31^2,35]. Therefore, 

T ^) = T W Xi0 (i) X20 (2)... X7v0(J v) 

By proposition 3.12 in [31] we get 

\\T' (N) 1 1 = | |TW x 1 0«x 2 0< 2 > • • • x N O W 1 1 = \ \T {N) 1 1 



Proposition 7 : If a multipartite pure state \ip) is subjected to a local measurement 
on the /cth qubit giving outcomes ik with probabilities Pi k and leaving residual iV-qubit 
pure state \4>i k ) then the expected entanglement ^2 ik Pi k Er(\(pi k )) of residual state is not 
greater then Er(\tp)). 
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5> fc £ r (|0 < E T (m. (26) 

Proof. Local measurements can be expressed as the tensor product matrix D = 
® ® • • • £g> D( N ^ on the expanded coherence vector T [36]. The expanded coherence 
vector T is the extended correlation tensor T (defined below) viewed as a vector in the 
real space of appropriate dimension. The extended correlation tensor T is defined by the 
equation 

1 3 

p = Yl T iii2-i N °ii ®Oi2®---®Pi N (27) 

where 0" ifc G {/, cr x , <7 y , (X 2 } is the i^th local Pauli operator on the fcth qubit (a = I) and 
the real coefficients % li2 ... iN are the components of the extended correlation tensor T. 
Eq. (2) and Eq.(27) are equivalent with T 000 -o = 1 , Thoo-o = s^, ^ li2 ... iM oo-o = 

^ 2 'X M} '--- and = T £li N i M*-" * 0. £(*>; fc = 1, 2, • • - n are 4 x 4 

matrices. Without losing generality, we can assume the local measurements to be POVMs, 
in which case = diag(l, D^) and the 3x3 matrix is contractive D^ T D^ < I 
[36]. The local POVMs acting on a iV-qubit state p corresponds to the map p i — ► -M(p) 
given by 



where are the linear, positive, trace preserving operators satisfying L^L^ = I 
and [L-^,L-^] = 0. The resulting correlation tensor of M(p) can be written as 

T {N) =T^x 1 D^x 2 D^---x N D^ 

where is 3 x 3 matrix and D^ T < I. 

Action of POVM on fcth qubit corresponds to the map •Mfc(|V')(V'l) = Yli Mi k pM} 

where M lk = I ® ■ ■ ■ L<V ® • • • / , £ ijk A'f^Sf = 1 and L Sf] = with th * resulting 
mixed state J2i k Pij0»fc I > where \4>i k ) is the iV-qubit pure state which results after the 
the outcome ik with probability pi k . 

The average entanglement of this state is 

(A °)I|-1 



Y^Pi k E r {\^i k )((t>i k \) = J2Pi k \\ r \ { £ k 

ik i k 



ik 



Y,Pik\\D {k) T {k) m\\- 1 



ik 
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where, by proposition 3.7 in [31], D^T^(\ip)) is the fcth matrix unfolding [24] of x k 

D^ k \ Therefore, from the definition of the Euclidean norm of a matrix, \\A\ \ = y/Tr(AA^) 
[37] we get 

E^^(i^)(^i) = E^[ Tr (^ (fe) %(iV'))^)(^)) J D (fc)T )] 1 - 1 

ik ik 



i^S°n-i = ^» 



-\4>) 

because D^ T D^ < J, and ^2 ik p ik = 1. We have also used the fact that Euclidean norm 
of a tensor equals that of any of its matrix unfoldings. ■ 

As an example, we consider the 4— qubit state [38] 

\i/>)abcd = 4=(|0000) + |0011) + |0101) + |0110) + |1010) + |nn)). (28) 

A POVM {Ax, A 2 } is performed on the subsystem A, which has the form A x = U 1 diag{a, [3}V 
and A 2 = U 2 diag{\/1 — a 2 , a/1 — (3 2 }V. Due to LU invariance of we need only con- 

sider the diagonal matrices in which the parameters are chosen to be a = 0.9 and (5 = 0.2. 
After the POVM, two outcomes |0i) = Ai\tp)/y/p^ and \(j> 2 ) = A 2 \ip) / 'y/p~2 are obtained, 
with the probabilities as p\ = 0.5533 and p 2 = 0.4467. We find 

E T (m = 0.7802, E T {\4>i)) = 0.0725/ Pl , E T (\<P 2 )) = 0.0436/p 2 . 

This gives, 

E T (m)-\p 1 ET(\<f> 1 ))+p 2 E T (\<f> 2 ))}= 0.6641 > 0. 

This is to be contrasted with the similar calculation in [38], with the same state \ip) in 
Eq.(28) and the same POVM given above. 

Proposition 8 : Let be an 7V-qubit pure state. Let p denote the reduced density 
matrix after tracing out one qubit from the state \ip). Then 

ll^ll < \\r^\\ 

with equality only when \ip) = \(f>) ® \x) where \x) is the state of the qubit which is traced 
out. 
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Proof, we prove this for a special case whose generalization is straightforward. Let 
\^) = a\b 1 ---b N ) + b\b' 1 ---b' N ); |a| 2 + |6| 2 = l. 

Here and are the eigenstates of operating on the ith qubit. Now consider 
set of TV-fold tensor products of qubit operators {a a }, a = 1,2,3, namely S = {a ai <g> 
<7 Q2 ® ■ ■ ■ <S> o- aN }, • • -a N = 1, 2, 3. Choosing a 1: • ■ ■ , a N = 3 we get a 3 <g> a 3 <g) • ■ ■ <g) 
o~ 3 \bi • • • &jv) = =b|6i • • • 6jv) ■ We can choose an operator from S, denoted B, such that 
B\b\ ■ ■ - b N ) = ±\b[ • • -fc^v). If -B contains q < N a x operators we can replace k < q of 
them by a y operators. We denote the resulting tensor product operator by B k (B = B). 
We have, B k \h ■■■b N ) = ±(i) k \b' 1 ■ ■ ■ b' N ). Then, 

(&i • • - b N \o- 3 (g) • • • <g> cr 3 |6i • • - b N ) = ±1 = (b[ ■ ■ -b' N \a 3 <g> ■ ■ ■ <g) a 3 \b[ ■ --b^) 

(b 1 ---b N \B\b' 1 ---b' N ) = ±l = (b' 1 ---b' N \B\b 1 ---b N ) 
(b' 1 ---b' N \B k \b 1 ---b N ) = ±(i) k 
(b 1 ---b N \B k \b' 1 ---b' N ) = ±(-i) k 

Now, 

t ai -a N = {ip\cr ai <g> • • • <g> <J aN 
= |a| 2 (6i • • -b N \a ai <g> • • • <g> <7 Qjv |&i • • • b N ) + |6| 2 (&i • • • b' N \a ai ■ ■ ■ ® a aN \b[ ■ ■ ■ b' N ) 
+a*b(b 1 ■ ■ -b N \a ai ® • • • <g) <7 Q J&i • • -6^) + afc*(&i • • •b' N \a ai ® • • • <g> cr ajv |&i • • -&at) 
The nonzero elements oft ai ... ajv are t 33 ... 3 = ±|a| 2 ±|5| 2 , t B = ±ab*±a*b = ±2\a\\b\cos(4> a — 

t Bk = ±(i) k ab* ± (-i) k a*b = I ±l\ a \ !J! C0S( ^ - t b \ * * 18 ™ 
^ w v ' \ ±2\a\ \b\sm((f) a - 4> h ) if k is odd 

We get Y^k=o (2I) e l emen ts with cos((p a — 4>b) an d Sfc=o (2*1+1) e l emerL ts with sin((p a — 4>b) ■ 
If q is odd (for the given state] V')) the number of cosines and the number of sines are 
equal. When q is even the number of cosines exceeds by 1. Finally we get 



|^ ) || 2 = (±|a| 2 ±|6| 2 ) 2 +4|a| 2 |6|W(0 a -0 b )^ [Ij+Mal^sin^-^ , 

fc=0 ^ ' k=0 ^ 



Note that, using |a| 2 + |6| 2 = 1, it is easy to see that ||7j^|| > 1, showing that Er > 0. 
Next we consider 

IV)^! = |a|>i---6 w )(6 1 ---6^| + |6| 2 |6;---6' w )(6;---6 / w |+a6*|6 1 ---^)(6;---^| 

+a*b\b[---b' N )(b l ---b N \ 
and trace out the Nth. qubit to get the N — 1 qubit reduced density matrix 

p = |a| 2 |&i • ••b N - 1 )(b 1 ■ • • 6jv-i I + \b\ 2 \b[ • ••b' N _ 1 )(b' 1 ■ •■b' N _ 1 \ 
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+ab*\b 1 • • -ftjv-iX&i • ■■b' N _ 1 \(b N \b' N ) + a*b\b[ • • • • •& JV _i|(&' JV |& JV ) 

Now 

*ai- ajv _i = Tr(pa ai <g> <j Q2 (8) • • -a^J = |a| 2 (&i • • -6jv-ikai <8> • • • <8> <T ajv _il&i • • 

+ I & | 2 ( & 1 • ' -6iV-lkai®- • •®<Ta w _ 1 |6i • • A-lK^^l • • -&JV-lkai®- • -^a^Jfti ' " ' V N -l) ( b N\b' N ) 

+ab*(b[ ■ ■•b' N _ 1 \a ai <g> • • • <g> <7 Qjv _il&i • • • &jv-i)(&jv|&jv) 

We have for iV— 1 tensor product operators o-3®o"3<8>- • • (H>cr 3 1 61 • • • &at_i) = ±|&i • • • 6jv-i)- 
We construct the operators £> and D k corresponding to B and B k acting on //-I qubits. 
We then get D|6i • • • &jv_i) = ±|&i • • •& , JV _ 1 ) and D fc |6i • • -6jv_i) = ±(i)*|&i • • Now > 
the nonzero elements oiTp N ~^ are t 33 ... 3 = ±|a| 2 ±|5| 2 , t D = ±ab* (b N \b' N ) ±a*b(b' N \bN) = 
2\a\\b\\{b' N \b N )\cos((f) a - 6 - a), 



fc w x 1 Nl v y at 1 \ ±2\a\\b\ \(b' N \b N )\sin((/) a - (f) b - a) if fc is odd 

Finally we get 

IIT^H 2 = (±|a| 2 ± |6| 2 ) 2 + 4|a| 2 |6| 2 |(^|^)|W(0 a - cf> b - a) £ (l 



k=0 



2k 



q 

2k + 1 



+4|a| 2 |&| 2 |(^M| 2 S m 2 (0 a - <t> b - a) £ 
where q' < q is the number of o\ operators in D. Since |(&jv|&at)| 2 < 1 we see that 

ii^ir^iiCn 2 

equality occurring when \b]y) = \b' N ) in which case \ip) = \(f>) <g> |6jv) ■ It is straightforward, 
but tedious to elevate is proof for the general case 

W) = ^2 a ai-a N \b ai • • -b aN ), OLi = 0, 1 
ai---a N 



Basically we have to keep track of B type of operators, where r is the number of 

AN 



terms in the expansion of in order to obtain all nonzero elements of When iVth 



particle is traced out, the corresponding elements of Tp N ^ get multiplied by the overlap 
amplitudes, which leads to the required result. ■ 

Continuity of E?: We show that for iV-qubit pure states ||(|^)(^| — |</>)(</>|)|| — > =>- 
Er{m - E T {\4>)) ^0 
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Proof. 10X01)11-0 



But \\T$> 
Therefore I 



\% 



2 \<t>) I 

(JV) 



> 



\%. 



£ r (|V» - £t(I0» 



0. 



l^> 



|0> M 
\ 2 W I 



^> 1 



_ iit^i 
- IM|0> I 



A. Entanglement of multiple copies of a given state 

LU invariance. We show that Er for multiple copies of iV-qubit pure state \ip) is LU 
invariant. Consider a system of N x k qubits in the state \x) = <S> <E> • • • <S> (k 
copies). It is straightforward to check that [24] 

t (jv) _ T {N) T {N) r {N) (C)Q) 

This implies, in a straightforward way, that 

ll i lx> II ~ W 2 W II ' 

Since by proposition 6 ||77^|| is LU invariant, so is II^a^II- 

Let \ip) be a iV-qubit pure state and \x) = <E> \ip)- Then Er(\x)) is expected to 
satisfy 

E T (\x))>E T (m 



We again use the fact that 



r {N) _ T (N) r {N) 

2 \x) ~ 2 W ° 2 W 



which gives 



i lx> II - I^W I 



Since H^ll >1 we get ll^ll > H^Hor, 

E T (\x))>E T (m. 
Superadditivity : We have to show, for iVqubit states \ip) and \<p) that 

Erm®\<t>))>E r m) + E T {\<t>)). (30) 
We already know that for \x) = \^>) <8> |0) 

iitWii - iit^h iitWii 

H J |x> II - II^W II H^> II 
Thus Eq. (30) gets transformed to 

\\T (N )\\ ||T (JV) II _ 1 > \\T (N) \\ -4- IIT (iV) ll _9 
IM W II II ^ IM W II + 11^) II ^ 

which is true for ||7j'^' ) || > 1 and ||7j'^' ) || > 1. ■ 
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B. Computational considerations 

Computation or experimental determination of Er involves 3^ elements of so 
that it increases exponentially with the number of qubits N. However, for many impor- 
tant classes of states, Er can be easily computed and increases only polynomially with 
N. We have already computed Er for the class of N qubit W states, GHZ states and 
their superpositions. We have also computed Er for an important physical system like 
ID Heisenberg antiferro magnet. For symmetric or antisymmetric states is super- 
symmetric, that is, the value of its elements are invariant under any permutation of its 
indices [24]. This reduces the problem to the computation of |(iV + l)(N + 2) distinct 
elements of T^ N \ which is quadratic in N. 



C. Entanglement dynamics : Grover algorithm 

We show that Er can quantify the evolution of entanglement. We consider Grover 's 
algorithm. The goal of Grover's algorithm is to convert the initial state of N qubits, say 
|0 • • • 0), to a state that has probability bounded above 1/2 of being in the state \ai ■ • • ajv), 
using 

U a \b 1 ---b N ) = (-l) us ^\b 1 ---b N ) 

fewest times possible. Grover showed that this can be done with 0(V2^) uses of U a by 
starting with the state 

"" -l 



x=0 

where 



H 



M 1 -i 



and then iterating the transformation H® N U a H® N U a on this state [23]. The initial state 
is a product state as is the target state, but intermediate states ip{k) are entangled for 
k > iterations. Figure 6 shows the development of Er(\ip(k))) with number of iterations 
k, for six qubits. The values of k for which Er vanishes are the iterations at which the 
probability of measuring \a± ■ ■ - ajv) is close to 1. Thus Er can be used to quantify the 
evolution of a iV-qubit entangled state. 

V. EXTENSION TO MIXED STATES 

The extension of Er to mixed states p can be made via the use of the convex roof or 
(hull) construction as was done for the entanglement of formation [16]. We define Er(p) 
as a minimum over all decompositions p = P« I V 5 *) (V'i I m t° pure states i.e. 

E r (p)= min V ft £ T (|^)). (31) 
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k 

figure 6. Entanglement in Grover's algorithm for six qubits as a function of number of 

iterations. 

The existence and uniqueness of the convex roof for E? is guaranteed because it is a 
continuous function on the set of pure states [39] . This entanglement measure is expected 
to satisfy conditions (a), (b) and (c) given in section IV and is expected to be (d) convex 
under discarding of information, i.e. 

Y^PiET(Pi)>Er(Z,PiPi)- (32) 

i i 

The criteria (a)-(d) above are considered to be the minimal set of requirements for 
any entanglement measure so that it is an entanglement monotone [29]. 

Evidently, criteria (a) and (b) are satisfied by Ej-(p) defined via convex roof as it is 
satisfied by E? for pure states. Condition (d) follows from the fact that every convex 
hull (roof) is a convex function [40]. We need to prove (c), which is summarized in the 
following proposition. 

Proposition 9: If a iV-qubit mixed state p is subjected to a local operation on % th 
qubit giving outcomes k with probabilities pk and leaving residual N qubit mixed state 
Pk, then the expected entanglement J2 k PkEr(Pk) of the residual state is not greater than 
the entanglement E r (p) of the original state. 

J^PkEriPk) < MP) 

k 

(If the operation is simply throwing away part of the system, then there will be only one 
value of k, with unit probability.) 
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The proof follows from the monotonicity of 3r(\ip)) for pure states that is propositions 
6, 7 and 8. Bennett et al. prove a version of proposition 9 in [34], which applies to any 
measure satisfying propositions 6,7 and 8. Thus the same proof applies to proposition 9, 
so we skip it. 

Note that any sequence of local operations comprises local operations drawn from the 
set of basic operations (i)-(iv) above, so that proposition 9 applies to any such sequence. 
Thus we can say that expected entanglement of a iV-qubit system, measured by Ej-(p), 
does not increase under local operations. 



VI. A RELATED ENTANGLEMENT MEASURE 

We consider following entanglement measure. Consider 

£ r (m = log 2 \\T^\\ = log 2 (E r (m + l) 

where is the iV-way correlation tensor occuring in the Bloch representation of 

Proofs of propositions 1 through 8 easily go through for £j-(\ip))- We prove continuity 
as follows. 

Continuity of £ r (\ip)) : We have to show, for two iV-qubit states and \<f>) that 

\mm - - o =► \e r m - mi</>»| - o 

we have, ||(|^| - \<t>){<t>\)\\ - =► \\T^ - T^\\ -> 

But||TW-T(f)||>|||rW||-||T(f)||| 
ii |V>> \<t>) ii — I ii \i>) ii n 10} 11 1 

Further, whenever |||7ft°|| > 1 and \T^\\ > 1 

we have \\\Tff\\ - \\7ff\\\ > \log 2 (\ | |) - log 2 {\ \1™ \ \) \ 



o 



s r (m-£r(\m 



0. 

However, £q-(\ifj)) has the added advantage that it is additive ( while Ej-(\ip)) is su- 
peradditive). Indeed, from section IV A we see that for k copies 

Similarly £ T {\^) <g> |0» = 5 T (|^» + 

The extension of £r(\ip)) to mixed states via convex roof construction is similar to that 
of Er(\ip))- Thus £r(\ip)) has all the properties of Er(\ip)), with an additional property 
that £r(\^>)) is additive, while Er(\ip}) is super additive. 



VII. CONCLUSION 
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In conclusion, we have developed an experimentally viable entanglement measure for N- 
qubit pure states, which passes almost all the tests for being a good entanglement measure. 
This is a global entanglement measure in the sense that it does not involve partitions or 
cuts of the system in its definition or calculation. This measure has quadratic compu- 
tational complexity for symmetric or antisymmetric states. Computational tractability 
is not a serious problem if N is not too large and the measure can be easily computed 
for systems comprising small number of qubits, which can have many important appli- 
cations such as teleportation of multiqubit states, quantum cryptography, dense coding, 
distributed evaluation of functions [41] etc. However, finding other classes of states for 
which Ef can be computed polynomially will be useful. It will be very interesting to seek 
applications of this measure to situations like quantum phase transitions [11], transfer 
of entanglement along spin chains [42], NOON states in quantum lithography [43] etc. 
Finally, we have extended our measure to the mixed states and established its various 
properties, in particular, its monotonicity. We may also note that neither its definition, 
nor its properties depend, in an essential way, on the fact that we are dealing with qubits, 
so that this measure can be defined and applied to a general iV-partite quantum system. 
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